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Abstract 

A periodic box-ball system (pBBS) is obtained by ultradiscretizing 
the periodic discrete Toda equation (pd Toda eq.). We show the relation 
between a Young diagram of the pBBS and a spectral curve of the pd 
Toda eq.. The formula for the fundamental cycle of the pBBS is obtained 
as a coUoraly. 

1 Preface 

A cellular automaton (CA) is a discrete dynamical system which consists of an 
array of a number of cells. Each cell allows for finitely many states which change 
into new states in discrete time. Usually the rule of time evolution with which 
the system is equipped is quite simple, and CA are often investigated as simple 
models for natural or social phenomena. The box-ball system (BBS) is one type 
of CA, represented by finitely many balls and countably many boxes arranged 
in a line. 

In this paper, we study a periodic box-ball system (pBBS), which is a BBS 
with a periodic boundary condition. The pBBS is closely related to integrable 
nonlinear equations. In fact, the pBBS has soliton-like solutions and a large 
number of conserved quantities. Moreover, the pBBS can be obtained from 
integrable equations by the method of 'ultradiscretization'. 

This relation gives us a new method to describe the behaviour of a pBBS. 
If the initial-value problem of integrable equations related to the pBBS is solv- 
able by some analytical method, the initial-value problem of pBBS itself is also 
solvable, as the solution of pBBS is obtained from the solution of the integrable 
equations by ultradiscretization. 

The present paper is organized as follows. In section[2l we introduce the def- 
inition of the pBBS and the pd Toda equation. These two objects are connected 
each other through 'ultradiscretization'. We define the conserved quantities of 
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these two systems and state a main theorem (Theorem I2.3P which yields direct 
relation between the spectral curve and the Young diagram. Section [3] is spend 
to prove theorem 12.31 In Section]!] and ]5] we give the solution of the initial 
value problem of pBBS and derive the fundamental period for it, as a corollary 
of theorem 12.31 

2 Periodic box-ball system and periodic discrete 
Toda equation 

2.1 pBBS 

Let us consider a one-dimensional array of L boxes. Let Q be the total number 
of balls, such that Q < L/2. Each of these boxes is either empty or is filled with 
a ball. Since we are interested in the periodic case, the A^-th box is adjacent to 
the first box. The time evolution of this system is: 

(i) In each filled box, create a copy of the ball. 

(ii) Move all copies once according to the following rules. 

(iii) Choose one of the copies and move it to the nearest empty box on the 
right of it. 

(iv) Choose one of the remaining copies and move it to the nearest empty box 
on the right of it. 

(v) Repeat the above procedure until all the copies have been moved. 

(vi) Delete all the original balls. 

It is not difficult to confirm that the resulting state does not depend on the choice 
of the copies. This dynamical system is called the periodic box-ball system, or 
pBBS. Figure JT] shows an example of the pBBS and its time evolution pattern. 
The last entry is considered to be adjacent to the first entry. The pBBS is usually 
regarded as a dynamical system of a finite sequence with periodic boundary 
condition. Let us denote an empty box by '0' and a filled box by '1'. 

Let be the number of groups of consecutive 'I's at t = 0. (Clearly, N 
is also the number of groups of of consecutive 'O's at t = 0). This number N 
does not change under the time evolution and it corresponds to the number of 
solutions in the pBBS. We introduce dependent variables Q*, Wj(j — 1, . . . ,N,t 
e N), as in figure ]2] 

At t — 0, choose one of the sets of consecutive 'I's and denote the number 
of 'I's by Qq. Next, looking to the right, denote the number of 'O's in the 
nearest set of consecutive 'O's by Wg . Then, looking to the right, denote the 
number of 'I's in the nearest set of consecutive 'I's by Ql- We continue to 
define Wi,Q2, ■ ■ ■ ,Q%,Wpf in a similar manner. Since our system has the 
periodic boundary condition, it follows Q% — Qq, Wp^ — WQ,...etc. In the 
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Figure 1: Time evolution rule for pBBS 
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Figure 2: The definition of Q*j and Wj 
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following, we always use the convention that the position j is defined in Zjv 
(i.e. Qj_|_Ar = Q], W/j+AT = Wj). At t = 1, to define Qg, Wg^, . . . ,etc., one needs 
to choose one of the sets of consecutive 'I's; The set of consecutive 'I's whose 
leftmost entry was updated from the '0' of the first set of consecutive 'O's will 
be called Qq. In general, Qp^^ is defined as the number of entries in the set of 
consecutive 'I's whose leftmost entry was updated from the '0' of the first set 
of consecutive 'O's at t. 

The following formulae describe the time evolution of the pBBS. 



Ql+' = min[iy/,X* + Q,*] (1) 



w!+' = Ql+i + w!-Q*^' (2) 



Xj = max 

fc=0,...,Af-l 



(3) 



where it follows that 



N N 



J2qI<Y.^! (4) 



i=l i=l 

due to the condition Q < L/2. 

The main feature we use to solve the initial value problem of pBBS is the 
correspondence between the pBBS and the periodic Toda equation. 

Definition 2.1 The periodic Toda equation (pd Toda equation) is given as 

= n+v^-v:^ (5) 
vr - (6) 

i 

with the boundary condition 

The following proposition shows the essential relation between the pBBS and 
the pd Toda equation. 

Proposition 2.1 ([Ij) Suppose that the pd Toda equation in definition\2A\ sat- 
isfies the condition 

N N 

o<n^/<n^*- (8) 

1=1 

// the pd Toda equation has a one parameter family of solutions Ijis) and 
Vj (e) (for parameter e ), then if 

g* = lim log /,*(£) and W*. = lim -elogV;*(£) 
exist, they satisfy the equations Q), (0], and 
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Proof. Substituting ([S]) to (O, we have 



Since if^l satisfies the same equation, 



jt+i 



jt i/t 



Repeating this procedure, we get the following equation due to the periodic 
boundary condition. 



/*+! = /* + V' 



n-i+yu 



TtTT- 



This is a quadratic equation of . The two solutions are 



rt+l 



and 



Tt+l 



= V' + If- 



1 



(9) 



(10) 



The first one does not satisfy the condition The other one gives the time 
evolution for /*. 

Now, we calculate the ultradiscrete limit of ^ and (jlOp. 



We put /* — exp 



is) 



Vl — exp 



and take a limit e +0. 



By virtue of the fact that 

lim -e log {e""/" + e~''^^) ^ min [a, b] 



and 



AT 



N 



1 - 



''l ^2 ■ ■ ■ 



Jt Tt Tt 



0, 



> lim e log 

i=l 1=1 

it is a straightforward result that JS]) yields ([2]), and ^ yields Q if Q\,Wl 
exist. ■ 



We shall use this proposition to solve the initial value problem of the pBBS. 
Our strategy can be summarized as follows: 
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(1) For given initial data Q°, (j = 1, 2, 
with the pd Toda equation as 



. , N), we associate initial values 



/" — exp 



K° = exp 



(11) 



(2) Then, we solve the initial value problem of the pd Toda equation by the 
inverse scattering method. The solution {lj{e),Vj{e)} depends on the 
parameter e. 

(3) In principle, by prop 12. 1[ the solution to the pBBS is obtained as 

Q*= hm -e log /*(£), W;= lim -elogl^*(£). (12) 



Remark 2.1 One does not always have to define /j* and Vj' as fll]). Indeed, 
one could define these numbers freely on the condition that 

g°= lim -£log/,°(e), W° ^ lim -£logV^°(e). (13) 



2.2 Solution of pd Toda equation 

The initial value problem of the pd Toda equation was solved by the algebro- 
geometric method [1]. In this paper, we omit most of the details of the method 
and give only the solution. 

Let C be a hyperelliptic curve of genus g, and define the base of i?i(C, Z) as 
in figure [3] We denote the normalized 1-form on C by {uJi} f^i, and the period 




Figure 3: Canonical basis of Hi{C,'Z). (Case for g — 3.) 
matrix of C hy B ^ ( J, tUj) 

Remark 2.2 {uJi} is normalized 44> tOj — Sij yi,j- 
A hyperelliptic curve C of degree g can be exprresed as 

= G(A), 
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where G{X) is a polynomial in A of degree 2g + 1 or 2g + 2. Any holomorphic 
differential on C can be rewritten as 

As-^dA Xa-^dX dX 

Cg-l \-Cg-2 1 h Co , Co, . . . , C„_l G C. 

M M M 

Let us define the complex constants Cj^kij = 1, 2, . . . , (7, k = 0,1, . . . , g — 1) as 

As-i dX A9-2 dX dX 

l^j = Cj,g-1 1- !-••• + Cjfi . 

M M A* 

The Abelian mapping and the theta function are the most important tools 
in the method. We define the quotient space /{Z^+BZ^) obtained by the 
equivalence relation 

For a fixed point Po £ C, the mapping 

a; = / toi,..., 

IPa \JPo JPo 

is a well-defined Abelian mapping. The Abelian mapping is usually denoted by 



C3P^ ( uj=[[ u;i,...,[ cOg ) e CV(ZS+BZS) 

JPo \JPn JPo / 



A{P)= f u>. 

JPo 

The Abelian mapping of a divisor D = J2i n-iPi is defined by the formula 

fPi 

A{D) = y^ni w. 

i 

Definition 2.2 Let B be a g x g matrix which satisfies the relation 

B = B* and ImB > 0. 

Then the theta function 9{z,B) for z gC^ is defined as the holomorphic func- 
tion 

6{z,B)= ^ exp(7rin*Bn + 27rin*z). 
Remark 2.3 The theta function 9{z, B) satisfies 

0{z + ej,B)=0{z,B), ej = {0,0,- ■■ ,0)*, 

and 

0{z + bj,B) = e-2^^^^-'^'^"^(z, B), bj = (Bij, ■■■ , Bg,y. 
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In our algebro-geometric method (i.e. inverse scattering method), we use 
these functions and Abelian integrals on some hypereUiptic curve C to describe 
the solution of pd Toda equaiton. 

To define the curve C associated to the initial condition {/°,V^°}^^^, we 
prepare two sequences {xn}nei and by (fTil) and ([T~ 



Xn+i = {X- {1^+, + V^)}X^ - m)x.n-l 



(14) 
(15) 



Let C be a hypereUiptic curve defined by 

Ai2 = A(A)2-4m2, (16) 

where 

N 

^{\)=XN+i+yN. m^ = \liy^. (17) 

1=1 

Note that A(A) is a monic polynomial in A of degree N . The genus of the 
hypereUiptic curve C is iV — 1 (=: g). 

Seperately form the definition of C, we define iV — 1 complex numbers /ij, 
(j = 1, 2, . . . , — 1) as the roots of 

VN+i = 0. (18) 

Note that un+i is a polynomial of degree iV — 1, the highest coefficient of which 

is -/oyo. 

Theorem 2.2 (^J) Let C he a hypereUiptic curve \16]) . and define the canon- 
ical basis of Hi{C,1j) as in figure\^ Then, the solution of the pd Toda equation 
0) & (0) is expressed as follows. 



3=0 ]=1 

d 9(nr + tu + c) 

^.^/^'^-^d^^°g g((n + l)r + ^:. + c) ' 

a 

where r = A{oo- - oo+), u = A(0 - oo+), c(0) = A(oo+ - ^ P^) - K. 

0O+ is the point at infinity in the upper sheat of C , and oo^ is the point at 
infinity in the lower sheat. is the point in the lower sheat with A(0) = 0. K is 
a Riemann constant of C and Aq, Ai, . . . , Ag are the roots of A(A) = 0. 

And Pj is a point on C, which satisfies X{Pj) = ^ij. 
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2.3 Young diagram 



In this section, we briefly review the correspondence between box-ball systems 
and Young diagrams. A Young diagram is a collection of boxes as shown in 
figure m We define a Young diagram associated to a state of pBBS. 



Pi P2 P3 Pi 



























L3 





Figure 4: The Young diagram associated to the state in figure [5] 

Let us consider the pBBS which has A^-solitons f Section 12. ip . When we 
regard the pBBS as a dynamical system of a finite sequence of 'O's and 'I's, we 
can introduce the following operation which we shall call '10-elimination'. 

(1) For a given state, connect all the '10' pairs in the sequence with arc lines. 

(2) Neglecting the 10 pairs which were connected in the first step, connect all 
the remaining 'lO's with arc lines. 

(3) Repeat the above procedure until all the 'I's are connected to 'O's. 

(4) Eliminate the 'lO's in a state, and obtain a new sequence. 

Figure [5] shows an example of 10-elimination. The mark '|' means '0-soliton' 




000001111000011011010001110000000 

00000111000 1|1|0 1 1 

Figure 5: An example of 10-elimination. A 3-soliton system with two 0-solitons 
is obtained from a 5-soliton system. 

which has no entry but has a position. A 0-soliton appears when we eliminate 
a soliton of length 1. We can perform this '10-elimination' repeatedly and 
transform any A^-soliton system into a {N — fc)-soliton system with k 0-solitons. 
Note that k is the number of the shortest solitons in the 7V-soliton state. Note 
also that the 0-solitons do not move under the time evolution rule. 
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Let pi be the number of 10 pairs in a state of the pBBS, connected with arc 
hnes in the first step of 10-elimination, (i.e. after one ehmination). Similarly, 
we denote by pj the number of 10 pairs in the j-th step of 10-elimination. Note 
that pi > p2 > ■ ■ ■ > pi, where / is the number of the last step. The most 
important aspect of these integers pj is the fact that the series {pi,p2, ■ ■ ■ tPi} 
are conserved quantities for the time evolution of the pBBS |3]. Using this series, 
we can associate a state of the pBBS with a Young diagram with pj boxes in 
the j-th column (j — 1,2, ... , I) (see figure S]). Then let us denote the distinct 
lengths of the rows by {ii, L2, . . . , Ls}. Note that Li > L2 > ■ ■ ■ > Lg. 

The following is a main theorem in this paper. 

Theorem 2.3 Let C : — A(A)^ — 4m^ be the hyperelliptic curve associated 
with the initial value problem of pd Toda equation defined by and 
let 



I 







id V,r := e- 



(see section \2.2\) . Then all of the roots of A(A)^ 
positive. 

Let us denote these by 



4m^ = are simple and 



< Ao < A+ < Ar < A+ < • • • < A- < A+ (g = TV - 1). 



Moreover, 



lim elogA - 



the length of the [N ~ j)-th row of the 

Young diagram associated with the state {Q", W^} 



A state of pBBS 



/O-exp-QO/e, 
K?-exp-M/0/£ 



pd Toda equation dS])-® 
with initial condition {/°, V^} 



Conserved quantity Conserved quantity 

Ultradiscretization 



Young diagram 



Theorem 12.31 



Curve C -.n^ = A(A)-^ - 4to^ 



Figure 6: The relation between the conserved quantities of two systems, the 
pBBS and the pd Toda equation. 

By virture of theorem 12.31 the ultradiscrete limits of the data which are de- 
scribed hy [j — Q,l, . . . , g) can be expressed in terms of the Young diagram. 
In fact, several essential data for the pd Toda equation — fundamental period 
etc. — can be expressed by A^s only. Hence, we can express the ultradiscrete 
limit of these data by the information of the associated Young diagram. 
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3 Proof of theorem 12.31 
3.1 Ultradiscrete limit 

For convenience, we first state several lemmas which we will use in the rest of 
the present paper. 

Definition 3.1 Let e be a positive number, and let f{e) be a continuous real- 
valued function ofe. Let us introduce the symble ' -> ' which means 

( 3S>0 s.t. < e < (5 ^ /(e) > 
i^ = lim,^o£log7(e) 

Remark 3.1 By definition, f ~t>F <^ log/ ^ F/e (e 0). 

Definition 3.2 The sign '~' stands for the relation 

f^g^ f^F, g^G and F ^G. 

Remark 3.2 The sign is not equivalent with the usual sign ' ~ ' i.e. 
f ^ g <^ 3m, M > s.t. < m < lime^o \ f{e)/g{e)\ < M < +oo. For ex- 
ample, if fie) = e~^, g{e) — e^^, then f^g but f g. Precisely speaking, 
under the condition f{e), g{e) >0 for <e<i I, f ^g / ~5 holds, but the 
inverse relation does not necessarily hold. 

Lemma 3.1 Let /i, /2, . . . , fn, g be continuous real-valued functions of e with 

g = fi + f2 + --- + fN.Lf 

g~oG, fi-oFi, fN^FN, 
then G = max {Fi, . . . , Fn}- 

Proof. By definition, for some number 5 > 0, /i(e) , . . . , /jv(£) ,g{£) are all 
positive if e e (0, 6). It can be assumed that Fi > F2 > ■ ■ ■ > Fn without loss 
of generality. For the largest numbers Fi = F2 — ■ ■ ■ = F„i (1 < to < N), one 
can rearrange the index if necessary and assume 

flhf2h---hfm, 

where /i ^ fc 4^ 3C > 0, limg^o \ ki£)/h{e)\ < C < +00. So, 5 ^ /i is obvious, 
therefore, noticing remark [3?2l and Def l3.2[ G = Fi is proved. ■ 

We have the following obvious lemma 

Lemma 3.2 Let f,g be continuous real-valued functions ofe. 

f^F, g^G => fg^F + G. 

Lemma 3.3 Let f,g,h be continuous real-valued functions ofe with /(e) = 
g{e) + h{e). Iff-oF, \g\ -t>G, \h\ -o H , andG^H, then F = mux [G, H]. 
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Proof. Without loss of generality, one can assume G > H. Let 5 be a positive 
number which admits 



< e < ^ 



/(£) > 
\h{e)\<Cs\9{e)\ 



for some Cs which depends on only 6 and Cs ^ ((5-^0). 

Thus, < / < |g| + < (1 + Cs) \g\ for e G (0, 5), &nA F < G holds. On 
the other hand, the inequality < j^l < |/| + \h\ gives 

1 < 77 + 77 < 77 + for £ G (0, d). 
Wl Ifl l5l 

As C5 for decreasing S, it follows that F > G. ■ 



Remark 3.3 If we omit the condition "G ^ H', the claim of lemma \3.3\ becomes 

F < max[G,H]. 

Lemma 3.4 Let f{X,e) be a polynomial in X with real coefficients of degree 
N +1: 

/(A,e) = A^+i - kN{e)X + kN-iie)X^-' -■■■ + (-if+'koie) 

where kj{e) > (j = 0, 1, . . . , A^) and kj ~t> Kj. Then, the roots of the equation 
f{X,e) = 0, Ao(e) < Ai(£) < • ■ • < Ag(e), satisfy 

Aat-oATat, XN-i~t>KN-i - Kn, ...,Xo~{>Kq-Ki. 

Proof. The fundamental relation between roots and coefficients gives 

fcjv = Ao + Ai H h Xn, 

fcjv-l = A0A1+ • ■ • + Ajv-iAat, 

ko — XqXi . . . Xn- 

Using lemma 13.11 and 13. 2[ it is easy to prove the lemma. ■ 

3.2 Ultradiscretization of the polynomial A(A) 

In this subsection, we define and calculate the key parameters associated to 
an initial state of the pBBS denoted by Uj, {j = 0,1, . . . , N — 1) and Pk, {k = 
0,1,..., TV — 2). In subsequent subsections, the ultradiscrete limit of the solution 
of pd Toda equation is expressed by Uj and Pk- 

Let C : — A(A)^ — Am? be the hyperelliptic curve defined by p^ . Note 
that A(A) is a monic polynomial in A of degree N{= g + 1). 

We use the following two propositions without the proof. 
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Proposition 3.5 The roots of 

A(A) = 

are all real and positive. 

Proof. The proof is given in [S] . ■ 

Proposition 3.6 The 2N roots of the equation 

A(A)2 - = (20) 
are real and positive. All of these roots are simple roots. 

Proof. The proof is given in [S] . ■ 

Definition 3.3 Let us denote A(A) by 

A(A) = A^+i - ug\9 + ug_i\3-^ + i-iy+^uo. 

We define the real numbers C/j (j = 0, 1, . . . , 5), as 



Uj :— — hm elosui, 



or equivalently, Uj -i> ~ Uj 



To define Pk, let us consider the polynomial yN+i{X) ( (fT8|) ). Note that 
j/Ar+i(A) is a polynomial of degree — 1 and the roots of yN+i{X) = are 
Hk (fc = 1,2, . . . ,5). 

Definition 3.4 Let us denote yi\f+i{X) as 

VN+i = -ly? {a» - t^g-iA^-i + «g-2A^-2 ^ . . . ^ (-l)9«o} 

(see hl8\) ). We define the real numbers (fc = 0, 1, . . . , 5 — 1), as 

Pk := - lim elogwfc, 

or equivalently, Ufe — 1> — Pk- 

To calculate Uj, we need to prepare a few notations. 
Let the set A{N) be 



A{N): 



{ai - 1, ai, . . . , Ofc - 1, Ok} 



fc= 1,2,. 

Vi, j, Oi ^ Oj ~ I A i ^ j ^ Oi ^ a 



, S.U 



An element of A{N) is a subset of Z/A^Z, which consists of pairs of two consec- 
utive numbers. (In Z/NZ, we regard N and 1 are consecutive numbers.) For 
iV > 3, the number of elements of A{N) is equal to F/v + where Fn is 

the N-th Fibonacci number. {Fn+2 = Fn+i + Fat, Fi = 1, F2 = 2). 
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Proposition 3.7 The polynomial A(A) = Xf^^i + y^r is expressed as 

iji-l.Ju...,jk~l,jk)eA{N) 

where {ji - IJi, ...,jk - 1, jfc} U {zi, . . . ,iN-2k} = Z/iVZ, and Yj = -IjVj, 
X, = A- (/,+!+ V-)- 

The proof of proposition l3.7l is elementary though shghtly involved. We therefore 
defer the proof to the Appendix. Defining J^+i + Vi and 6j := liVi, we 

find: 

Proposition 3.8 The coefficients of the polynomial 

A(A) = Af+i - UgXS + Ug^iXs-^ + i-iy+^uo 

satisfy 

uo~> -{Qi + Q2 + --- + Qn) (21) 

and 

ug^ -mm[Q„W,]. (22) 
Equivalently, Uq = Qi + Q2 + ■ ■ ■ + Qn and Ug = min [Qi, Wi] . 
Proof. It is sufficient to prove 

Uo = i-bji)---i-bjk)ah---a.,^_^^ (23) 

(ji-iji,....]k-i,jk)eA{N) 

= hh.-.lN + V1V2...VN (24) 

and 

Ug = Il+I2^ V In + Vl + V2^ Wn. (25) 

is a direct consequence of proposition [X71 It remains to prove Substi- 
tuting Oi = + Vi and 6^ = /iVi to (|23p . two types of terms appear, namely 
those that contain Vklk for some /c (type (i)), and those do not (type (ii)). 

Among all the terms in (I^S)) . a contribution 'Vfe/fc' must come from the term 
which contains —hk or Ok-iOk only. For any term which contains —bk, there 
exists one term in which —bk is replaced by ak~iak in (j23p . Hence we can 
conclude that the summation of all terms of type (i) will cancel out. The only 
terms of type (ii) are /1/2 ■ ■ ■ In and V1V2 ■ ■ ■ Vn because the term of this type 
must come from 0102 ... a^r. ■ 

The ultradiscrete limit of ui, M2, . . . , Wg-i are also obtained in a similar man- 
ner. Let 

X = {A, I ^2,-1 - Qi, A2i^Wi, (/ = 1, 2, . . . , N)} 

and 



1 < cr(l) < cr(2) < • • ■ < cr(iV ~ k) < 2N 
(T{i) + 1 <a{i + 1), \fi 
cr(l) = 1 =^ (t(7V -k)^2N 
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Proposition 3.9 It follows that 

Uk -l> - min < ^ Acr(j) > (= -Uk). 

[{A„(.)}Gi3(fc,A') J 

Proof. From the proof of proposition 13.81 Uk can be obtained in the following 
way. First, calculate 

0,^0,2 ...a,„„,. (26) 

{n,Z2,...,iN-fc}CZ/AfZ 

And pick up the terms which contain no Vilis. Since ai = + Vi, a term that 
contains Vili+i cannot exist in ([26]) . Conversely, a term of length N — k which 
has neither Vili nor Vj/j+i necessarily appears in 



We can calculate Pq, Pi, . . . , Pg-i analogously. 
Proposition 3.10 The polynomial ?/jv+i(A) is of the form 

Oi-iji,...jfc-ijfc)e^'(7V) 

where X and Y are given in proposition\3.7\ and A'{N) = A{N) n 2(^/^^"'t^^) 
is a subset of A{N) of which element (ji — . . . ,jk — l,ifc) does not contain 
the number '1' G 1/NZ. 

We will prove this proposition in the Appendix. 

In the same way as in proposition 13.81 and proposition 13.91 we obtain: 

Proposition 3.11 It follows that 

Vk--^ - min < Y ^<t(») ( (= 

w/iere6'(A;,iV) =S(A;,iV)n2(^-{'32,H'i}) is a subset of B{k, N) which does not 
contain Q2 nor Wi . 



Proof. As in the proof of proposition 13. 9[ Vk can be expressed as follows 



Vk 



_{ii,i2,...,ijv-i-fc}C(Z/JVZ-{l}) 

Since oi ^ I2 + Vi, we obtain the proposition. ■ 

The following lemma can be obtained from proposition 13 . 91 by combinatorial 
arguments, which we will give in the Appendix. Let Zn be a set of an 7V-soliton 
state of pBBS. 
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]c/4 



=1 



U: 



2=0 



C/o=12 



Figure 7: Example of the interpretation of the [/(,s of lemma 13.121 

Lemma 3.12 Let x S he an N-soliton state, and let Uq, Ui, . . . ,Ug be the 

positive integers defined as in Def \3.3[ Then Uk is equal to the number of boxes 
below the k-th row in the Young diagram. 

Corollary 3.13 Let x G Zjq he an N-soliton state. We can obtain the integers 
Pq, Pi, . . . , Pg-i in Def \S.4\ by the following procedure. 

(1) Let Q^{x), Wj{x) {j — 0,1, . . . , g) be the integers defined in section [Kl\ for 
a state x G -Ejy. Consider another N-soliton state y G , with 

Q%y)^Q%x) for 3^2, 
WUy) = Wf{x) for i^l, 

Ql{y) » I, and W^{y) » 1. 

(2) Let pi be the length of l-th row of the Young diagram corresponding to 
y G Zn. Then, 

3-1 

P. = E 

1=3+2 



Pn=S 



Pi=5 



Figure 8: Example of a Young diagram corresponding to a state y. 
Proof. From proposition 13.91 and 13. IH it is obvious that 

lim iti+i = Vj. 
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Then, 



P, = — lim eloffw, = — lim lim = lim J/i+i. 



The fact /2, Fi ^ +0 <^ Q2, VFi ^ cxd completes the proof. 



3.3 Ultradiscretization of the curve /i^ = A(A)^ — Arm? 



To complete the proof of theorem 12.31 we consider the asymptotic behaviour of 
the Riemann surface C : fi^ = A(A)^ — 4™^ when e ^ 0. From p?|) . we easily 
obtain that — e~~ , or equivalently 

m-^-L/2, (27) 

where L is the number of boxes in the pBBS (see (fT7|) ). 

Recall that we have denoted the roots of A(A)^ — 4m^ = by 

< A(7 < A[!" < Aj" < Aj'" < • • • < Ag < A+. 

and the roots of A (A) =0 by 

< Ao < Ai < • • • < Ag. 

Note that A^^ < Xj < A+ {j = 0,1, . . . , g) . 

It is not easy to calculate the asymptotic behaviour of an Abelian integral 
on a general Riemann surface. However, as we shall see below, the problem 
of describing the asymptotic behaviour of C can be reduced to that of the 
degenerate case 

Co:a*' = A(A)2. 
By definition, it follows that 

A(A)2 - 4m2 = (A - A„ )(A - A+) ■ • • (A - Ag-)(A - A+), (28) 

A(A) = (A- Ao)---(A-Ag). (29) 
Clearly, equation ([28|) can be decomposed 

a 

A{\) + 2m = Y[{\-Xf^) (30) 

3=0 



and 



A(A)-2m = n(A-A7-(-')) (31) 
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where erf?) = I ^ ['1 ^ ^' a '^ \ and —(7(7) denotes the opposite 

^ ' \ - (j =5,5-2,5-4,...) 

sign to <7{j). 

By (HID, (gni) and dSD), we have 



A^. - a;« = -'"^ , (32) 



Lemma 3.14 It follows that 



Proof. Let uq be the constant term of the polynomial A (A) (see Def l3.3p . Since 
Qi + Q2 + ■ ■ ■ + Qn < L/2 (see section [^TT|) . we obtain 

uo ^ (uo ± 2m), 

from lemma [231 proposition 13. 8[ and (|27p . The proof then follows from 
(P^ . and lemma I 



As a corollary of lemma [3^ lemma [3.121 and lemma [3 .141 we obtain theorem 

[231 

The following lemmas and propositions are used to calculate the ultradiscrete 
limit which is expressed as the difference of XjS. 

Proposition 3.15 Under the condition 

j^k ^ Uj+i-Uj ^Uk+i~Uk, (34) 

the ultradiscrete limit of | Xj — | satisfies 

\Xj — A^l -i>J7j+i — Uj (j > /c), and 
I Aj - Afe I -i> Uk+i - Uk U < k), 

where Uj (j = 0, 1, . . . , 5) is the real number defined by Def \3.Sl 

Proof. From lemma 13.41 and lemma I3.14[ it follows that 

Xj ~t>Uj+i — Uj, and X"^ -t>Uk+i — Uk- (35) 

The assertion is proved immediately by virture of lemma [X5I and ([M[) . ■ 

Remark 3.4 The claim of proposition \3.15\ is also true without the condi- 
tion ^34^ . We can prove this assertion by using the fact that C/j+i — Uj {j = 
0, 1, • • • ,5) can be perturbed independently over the real numbers, as we can nat- 
urally extend the domain of the initial condition of the pBBS {Qn, W^} C N C 
]R>o. Continuity of the time evolution rule of pBBS justifies the 

argument using such small perturbations. 
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4 Ultradiscretization of the solution of pd Toda 
equation 

4.1 ultradiscretization of Abelian integrals 

It is not easy to describe the normahzed holoniorphic differential of a Riemann 
surface in general case. However, it is easy in the case of the degenerate curve 
Cq. In fact, the normalized holomorphic differential of Co is expressed as 



Theorem 4.1 It follows that 

' (e^O). 



bi 



To prove theorem 14. 11 we prepare several lemmas. Let be holomorphic 

differentials on C defined by 



-,o 



1 A(A) r 1 1 



d\. 



^ ' 27ri ^A(A)2 - 4m2 I A - A^ A - Aq 
In the first place, let us prove: 
Lemma 4.2 It follows that 

Jbi Jbi 

Proof. Let Xi, (i — 0, 1, . . . , 5 — 1) be real numbers which satisfy 



for example Xi := J Xl X.^^-^. Then we obtain 



(36) 



A+<X,<A,:;„ { ; (37) 



[' , [/. ' 1" (38) 

7a+ ^/A{X)^ - 4m2 I A - A, A - Aq J 
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Under the condition < A < Xi 



Uk^^ (A-A-)(A-A+; 




where 



(43) 



(44) 



(The existence of positive numbers I^^I^. in (j44p is proved from proposition 



13. 8[ proposition 13.91 and \21\. and the fact 

Q[! + Q? + • • ■ + 

For example, we can show 



2 



Aj, — Afc 



Ai — Afc 



< 



-t> 



2m 



AnAi ■ ■ ■ An 



, (e « 1) 



{Ql + Qi 



< 0, 



which assures the existence of /j, > 0.) 

Then there exist a positive number B' > such that 

1 + Oje-B'/'^) 

27ri 

l + 0(e'^'/") 
27ri 

-2 log 



^' i7(A-Ar)(A-A+) ^(A-Ao)(A-Ao+) 
2 log 



dX 



XT - a; 



A+ - a; 



I Xj- X^ ^ ^ / X, - Aq 

-^0 ^ ^0 





[A 




+ 2 log 









1 , -^i — A^ — Aq 

— log -, 

27ri X, - Ao A+ - A. 



Xi 



(45) 
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In the case of i ^ j, it is easy to show 

In a similar manner, it follows that 

f-Ar 



Xi 



OJ" 



(46) 



(47) 



Xi 



Equation (gS]), (gS]), and (gT]) complete the proof. ■ 
Next, we prove the following lemma, which completes the proof of theorem 

Lemma 4.3 It follows that 



Proof. Since ujj is a holomorphic differential on C, the Ricmann bilinear equa- 
tion [5] gives 



where 



J ai J hi J ai 

Then we obtain 

9 



If * ^ j, 



A" 



2 

27ri 



Jx- \ \/A(A)2 - 4to2 ^A(A)2 -4m2 J 



(48) 



(49) 



< 



1 + 0(e-^"/-') 



(A - AOdA 



K ^/(A-A-)(A-A+)(A-A-)(A-A+) 



(A - A,)dA 



7(A-Ar)(A-A+)(A-Ao)(A-A+) 



A+-A: 



|A. -Afc| + |A, -Aol 



(50) 

(51) 
(52) 
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where 

B" >0, and B'" > (53) 
The relation ([55)1 can be shown in a way similar to In the case oi i = j, 



1, 



>d\ (54) 
(55) 



TTiA- \ V^(A)^ -4m2 ^A(A)2 - ^w? 

'^r ./(A-Ar)(A-A+) 



where F > 0. Substituting ([52]) and ([551) to then 



By theorem l4.1[ we can calculate the asymptotic behaviour in the limit e — > 
of the period matrix B = {Bij)Ki j<g associated with the hyperelliptic curve 
C : = A(A)2 _ 4^2. 



Bj 



1 



1 



1 

7ri 



1 



27ri I A — A; A — Aq 



1 A,; — Aj Art — Ao 



dX 



\ - Ao Aq - Aj J 
Other parameters can be calculated similarly: 



1 



1 



1 



-oo 27ri [A - Xj 



A -An 



dX 



1 , A, 
= log IT' 

ZTTl Ao 



An — Ao 



Ao — Ao 



^ 1 1 Aq — Aj 



(56) 
(57) 

(58) 
(59) 

(60) 
(61) 



E 



1=1" ^0 



E 

i=l 



.° + k, 
1 i-1 



1=1 -""i 



(62) 

(63) 
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Using the formula for the Ricmann constant corresponding to the hyperelliptic 
curve 



([63]) becomes 

a 



E 



E 

i— 
9 

E 

1=1 

ff 

E 



(64) 



1 

'27ri 



A — \j A — Aq 



i=l 



1 - Aj - Aq 1 

log ■ 1 

27ri A,- - Ao Mi - 2 



(65) 
(66) 



Now, using remark im proposition 13. 15[ and (|32)). 



R . ~ — log 



4m^ 



TTi (Aj - Ao)2 nfc>l(Afc - Ao) Uo<k<ji^J - Afc) nj</c<g(A/c - Aj) 



TTi e ^ — ' ^ — ' 

j<k<g 



k>l 



-- ( 2Af - (j + l)C/,+i + (j + 2)C/, + C/i 

TTI £ 



and for i > j 



B,. 



2m{\i — Xj) 



— log ■ 

TTi {Xj - Ao)(Ai - Ao) nfc>i (Afe " Ao) 

1 1 
— -I 

TTI e 

1 1 

7ri e 



:(M - (t/,+1 - U,) - J2 (Uk+i - Uk)) 
fe>l 

1-L(m -Uj+i + Uj + Ui 

me \ / 



(68) 
(69) 

(70) 



where m-i>M(= -L/2), and C/g+i = 0. 
Similarly, 

1 1 



{r),^-~[M-U,+, + U, + U, 



(71) 
(72) 
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and the j-th elements of Cq 



0+ 9 

j^l -'mo 



Coj 



E 

9 

E 

4=1 



1 1^ (A^"A,)(Ar-Ao)(M.-A,-) 1 
^(A+-Ao)(A--A,)(m.-Ao) 2 



^ 1 , (A,--A+)(Ar-Ao)(^,-A,-) 
2^i ^Ao+-Ao)(A--A,)(m.-Ao) 



(73) 
(74) 



where we have chosen the branch log(— 1) = iri. Using the fact that fij [j = 
l,2,...,g) are the roots of yM+i{^) = 0, we can immediately calculate the 
ultradiscrete limit of all the terms in (|74p except (fii — Xj). (Note that -oPi — 
Po, fik-oPk - Pk-i, 



. , /ig -o - Pg-i.) Indeed, 

9 



n(A--Ao)^5](c/.+i-;7,; 



-Ui, 



n(Ar-A,) 



and 



n (A,-A,) X (A, -A-) X n(^'-^^) 



l<i<] 



j<i<g 



n(^^--^')X n lA^'-AI ^ X{(>^^-h) 
l<i<i -ll'^i I ' J I i<i<q 



-t> M-iUi-Uo), 



n (m. - Ao) 



i=l 



i=l 
l<i<9 

= -Po. 



Unfortunately, it is not easy to calculate the terms {^li — Xj). For the 

i=l 

time being, we treat these terms formally as 

9 
i=l 

We will prove a concrete expression for Sj in the Appendix. Thus, we obtain 
coj ^ ^ ( 9{M - C/,+1 + [/, + C/o) - Po - S,) . (75) 
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The fundamental decomposition of B is given as 

/Ai ... \ 





1 1 
2^ie 



A2 
V 



where 



and 



k + 
k 


\ M 


- (fc + 


l)Uk + kUk+i), 


(fc 


= 1,2,. 


■■l9l Ug+l = 




( 1 





... 0\ 




/I 





... 0\ 




-1/2 


1 


... 







1 


... 




-1/2 


1 


... 







-1/3 


1 ... 




1-1/2 





... l) 




u 


-1/3 


... l) 



Thus, using (|7ip . we obtain 



where 



^2 
?3 



i?-V=-(^i,^2,.. 
Ibi , 1 62 , 1 b. 






1 . . 





.. 


Vo 


.. 






5 . . 


. . -I- 











-1/.9 1/ 



2 ci 2 • 3 C2 3 • 4 C3 
I62 1 63 
3 C2 3 • 4 C3 

_163 
4C3 



5(5 ^ 




1 








1 










1 bg 


5 


+ lCg 



with 

6fe = C/o - (fc + l)C^fc + fcC^fc+i, Cfe = -M - (fc + l)Uk + kUk+i. 
On the other hand, ((67| . ([70| . and ((72| yield the important relation 

1 



-(bi 



5 + 1 
— (bi + bo 



where the bj is j-th column vector of the period matrix B. 
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4.2 Ultradiscretization of the theta function solution to 
the pd Toda 

In this subsection, we will calculate the ultradiscrete limit of the meromorphic 
function of the form 

^^•("^ ■■-^^''^ eiz-Ub[T\b,lBy -here. = (z,...,z,y (80) 

rather than the theta function itself, because we want to ultradiscretize ([1 
with ((79)) . We introduce the real matrix B° and the real vector z° hy B 
iB°, z = iB°z°. Starting from DefESl 



0{z,B) = exp (7rin*Bn + 27rin*2) (81) 

= exp(-7rn*B°n-27rn*B°2°) (82) 

= J2 expi-Tr{n + z°YB°{n + z°))expi7rz°*B°z°), (83) 



and 



d 

—0(z,B) = 27ri Y exp (7rin*Sn + 27rin*2:) (84) 



27ri exp (-7r(n + z°YB°{n + z°)) exp (7r;z°*S°z°05) 



Using these formulae, (|80p becomes 

6i,(2;°)6i(;z°-e)-6l(2:°)6i,,(z°-e) , , 

^^ T.n,m jnj -rrij) exp (-g(n +z°)) exp (-^(m + 2:°- e)) 
En,m exp (-ii'(n + z°)) exp (-7f(m + z° - e)) 

where e = (1/iV, . . . ,l/Ny, and i/(a;) = x*{ttB°)x. 

Since _B 0{e^^), we define r(a;) := lim^^+o (^c)- In (|87p . it turns out 
that the ultradiscrete behaviour of ^'j(z) is strongly dependent on the term 
(jij—nij). Recalling the fact that the period matrix must satisfy Im_B > 0, we 
find that r(a;) > 0, Va; G M^. Since F is a quadratic form over R^, we can order 
all the elements of x as 

< T{z° + n(i)) + Tiz° + m(i) - e) < r{z° + n^^)) + T{z° + m^^) - e) 

< T{z° + n(-^') + Tiz° + m(3) - e) < . . . , 

((nW, m,W) e Zf X Zs, i = 1,2,. ..)• 
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Let nf' and rrSf' be fc-th element of n^*) and m^^\ Then, the asymptotic 
behaviour of j{z) is described as 



■^j{z) ~ 27ri^ (nf - ) exp — (r(n(^) + z°) + r(m(*) + 2° - e) 

2—1 



Let 



z — e 



G^{z) r(n^*^ + + r(m^*^ + - e) - r(n^^^ + z°) - T{m^'^ + z° ~ e). 

Since F is a positive definite quadratic form over W ^ the set {a; e | |a;| = 
R} C W is bounded for any i? > 0. Thus, the set 



[i e N|r(nW 



r(m 



R} 



is a finite set for any R. We arrange all the elements of {Gi{z)}i^^ as 
< Gi(2;) = • • • = G^^l){z) < G^(i)+i(2;) = • • • = G^(2)(2;) < . . . . 
The relation (l88t becomes 



where /p = Eri!i\p_i)+i " "^j'O- 
Let 

q{j) := niin {p G N | /p 7^ 0}. 

To calculate the ultradiscrete limit of p^ . we recall the calculations in 
section Em and notice the following relation, 



00 crip) 

27ri^ ^ (nj'' — m^-'') exp 

p=l ;=a(p-l) + l 
00 

p=i i=<j(p-i)+i 

r, ■ r r G'o-(p)(2;) 

27ri > /pexp 
,(0 



Giiz) 



— m ^ j exp 



Gz(2) 



(90) 
(91) 



which is obtained from ujj ^ lo^ = 



coefficient c, g-i defined in Secl2]is found to be: 



27ri I A — Aj A — Ao 



d\. Hence, the 



27ri 



27 



Substituting these relations in (fTO|) . we obtain 

+ = ^(/° + l^°)-£(A.-Ao)-£^^*^(^) (92) 

i=0 i=0 j=l 

= Ea.-E(A.-Ao)-^^^v,^(.) (93) 

i=0 i=0 j = l 

= (g + l)Ao-^^^*,(^), (94) 
i=i 

where z — nm + tv + c(0). 

The following formula gives an answer to the initial value problem of pBBS. 



Theorem 4.4 The ultradiscretization of is given by 

min [Q*,+2, W^+,] = niin [u^ - [/i, {U, - U,+,) + G,(,)(^)] , (95) 

where Gq(^j){z) lim Gq[j){z) with z — nm + tu + c(0). 

Proof. By lemma 13. 3[ is obvious when Uq — Ui, and {Uj — Uj+i) + 

Gqi^j-){z), [j — 1,2,...,(7) are all distinct. In the general case, we have only 
to consider small perturbations as in remark Since we can make Uq — Ui 
and {Uj — Uj+i) + Gq(^j) {z) all distinct by perturbing Uj — Uj+i independently, 
we can conclude that ([95]) holds in the generl case by continuity of both sides 
of the equation. ■ 

Remark 4.1 Note that j{nm + ii/ + c(0)) does not change under the trans- 
lation 

tv 1-^ tu + bk, Vfc 

or equivalently 

t(B-%v^t{B-%u + l, 
where {B~^)k is the k-th column of B^^. 

5 Fundamental cycles of the periodic box-ball 
systems 

5.1 relative period 

Let Zjv be a set of A^-soliton states with no 0-soliton. We treat separately the 
set of A^-soliton states with 0-solitons, which is denoted by Z'^. A state x G Zjy 
is expressed as 

x = xiX2---Xl for Xi£{0,l}. 
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Using the translation map 

S : Zjsi Zjsi X1X2 ■ ■ ■ XL ^ X2 ■ ■ ■ xlXi 

which sends the first letter to the last, we define the set T{x) C Zn (for x G 
Zn), 

T{x) -.^{ye ZN\3meZ>o s.t. S'"(x)=y}. 

Let us denote the 10-elimination by El : Z^ ^ Z^q U lJn<iv -^n^ time 
evolution in the pBBS by T : Zn Zn ■ We also define V : Zn U Z'^ — * Zn 
as the map which acts as the identity on Zn and eliminates the 0-solitons in 

Remark 5.1 El is bijective. o El is surjective, but not injective. 

Definition 5.1 The fundamental cycle f{x) of x G Z^lor G Z^) is the min- 
imum positive integer p that satisfies T'p{x) — x. The relative period of x & 
Zjq{or G -Zj^), r{x), is the minimum positive integer q for which T'^{x) G T{x). 

Remark 5.2 r{x) \ f{x) for any x G Zpj because x G T{x). 

Remark 5.3 S and T commute. And S and El also commute. Consequently, 

T{x)^T{y)^T{¥.\{x))=Tmy)). 

Since the method we presented in section [2] can only be used to calculate the 
relative period of pBBS systems, the following claims are important. 

Lemma 5.1 Let x G Zjq be an N-soliton state of the pBBS. It holds that 
r(EloT"(a;)) T(r" oEl(x)), neE 

or 

r (Er 1 o T" (x) ) = r (T" o Er ^ (x) ) . 

Proof. Let Ql^ and [n = 1,2, . . . , N, t G N) be numbers defined by an 
A^-soliton state x (see section [^?T1 figure El) . Note that the equations (HJ-lHl) give 
Ql^^ and W^+^ from and W^. Obviously, when we replace Ql^ and W,* to 
Ql^ — 1 and — I, then these equations give Ql^^ — 1 and — 1, which 

means 

T(Elor(a;)) = T(ToEl(a;)). 

Using this formula, the former assertion is easily proved by induction. To prove 
the latter assertion, we start from the formar assertion. 

r(El o T" (a;) )= T(T" oEl(a;)) ^ 3m > s.t. El o r"(a;) S™ o T" o El(a;) 

^ Elor"oEr^(y) = S"oT"(y) 

^ r"oEri(y) = ErioS™or"(y) (96) 

where y = El(a;). Recalling remark [5T51 (ITO]) is equivalent to T" o EP^(y) = 
S™ o EP^ o T"(y), which completes the proof. ■ 
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Proposition 5.2 // a state x G Zf^ satisfies the condition El(a;) G Z'^_-^, then 
r(x) =/(El(a;)). 

Proof. Let x = El(x). Since El(a;) G x has exactly one 0-soliton. If 

TP{S:) = i, 

r{x) = T{Er\i)) r(Eri o TP(i)) ^ t{tp o Eri(i)) = r(rp(a;)). 

So, r(a;) < /(i). Conversely, if r(x) = T{T''{x)), it follows that 

r(i) r(El(a;)) = r(El o T«(x)) = r(T9 o El(2;)) = T{T'^{S:)). 

Since a 0-soliton does not move under the time evolution, the fact that x has 
exactly one 0-soliton leads to i = T*(i). So, r{x) > f{x). ■ 

Remark 5.4 In the vroof of vrovosition VST^ . we conclude thatT{x) = T{T'^(x)) 
x = T'^ix). This claim Jails in the case where there are more than two 0-solitons 
in the sequence o/El(a;), where are arranged symmetrically. In this situation, 

r{x) % f{i). 



call this symmetry 'internal symmetry'. Internal symmetry makes the prob- 
more complicated. We do not concider this symmetry in the present paper. 



We 

lem more 

The statement of proposition 15.21 can be generalized as follows. 

Corollary 5.3 // a state x G 2jv satisfies the condition El(a;) G Z*, where 
n < N, and x is without internal symmetry, then r{x) = /(El(a;)). 



By virtue of Cor l5.3l the fundamental period of the pBBS can be obtained from 
the relative period of the corresponding pBBS. 

5.2 Formula for the fundamental period 

Recall the definition of the Young diagram associated with the state of pBBS 
(Section [Q]) . We also define 

ni :— {the number of the rows of which length is L;}, 

I 3 

lo := L -'^2pj, Ij := Lj - Lj+i, N-j := Iq +'^2ni{Li - Lj+i). 

3=1 1=1 

The fundamental cycle of x G Z^ can be described by using the data of 
the corresponding Young diagram. In fact, the following formula gives the 
fundamental cycle of the pBBS system. 

Theorem 5.4 Let x G Z^ be a N-soliton pBBS without internal symmetry. 
Then 

fix) = LCM (Efl^fzl ^1^0 
V Islo Is-ilo hlf) 
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Remark 5.5 Recalling vrovosition \5.3l the formula in proposition \K^ is equiv- 
alent to 

r{x) = LCM — , . . . , -r-^, 1 

because one can obtain the Young diagram corresponding to Y o El(a;) by elimi- 
nating the first column of the Young diagram corresponding to x. 

Though this formula was first obtained by elementary combinatorial methods 
[4], we can obtain the same formula using a different method relying on the 
results of the previous sections. 

From (|19p . ([77)) . and remark 33 the relative period r{x) satisfies 

r(x)= LCM f-, -,lY (97) 



where are numbers defined by l[77|) . We use the following lemmas to prove 
theorem 15.41 

This lemma claims that there is a simple relation between the ultradiscrete 
limit of Riemann surfaces and Young diagrams. The proof of this lemma is 
given in the Appendix. 

Corollary 5.5 <;fe — c^fe-f i <^ the lengths of k-th and [k + l)-th rows from the 
bottom in the Young diagram are equal. 



Proof. 



1 &fe , 1 bk+i 1 bk+i 



fc + lcfc (fc + l)(fc + 2)cfe+i k + 2ck+i 



h bk+i 



Ck Ck+l 

Uo-(k + l)Uk + kUk+i _Uo-{k + 2)Uk+i + (k + l)Uk+2 



-M - (fc + l)Uk + kUk+i -M -{k + 2)Uk+i + ik + l)Uk+2 
<^ Uk — Uk+i = Uk+1 — Uk+2- 

Lemma [3.121 completes the proof. ■ 

The following lemma is almost trivial. We omit the proof. 

Lemma 5.6 Let p,p' ,p" , q,q' ,q" be integers, and p,p' ,p" are relatively prime 
toq,q',q" respectively. If 

9 _ 9^ ^ 9^ 
p p p 

then LCM {q, q') = LCM (q, q"). 

Proof of theorem 15.41 First, we prove the theorem for the case where 

i^3^'i^^^3■ (98) 
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Starting from 

r{x) = LCM (—,...,— ,1 ) (99) 



a 



2 2 2 2 \ 

= LCM , , . . . , , . . . , , 1 100) 

\'?i + Z^j=i^j ^1 — <r2 <r/c — <r/c+i '^a-^^'^a I 

by lemma 15.61 From ([77)) , 

I bk , ( 1 1 \ 



fc+lcfc V(fc + l)(fc + 2) k + 2J Ck+i 

{-M -U^){Uk+2-2Uk+i + Uk) 

(-M - [k + l)Uk + kUk+i){-M -{k + 2)Uk+i + (k + l)Uk+2) ■ 

On the other hand, by definition of Ij, Nj and lemma [3.12[ we obtain 

I, = - L, = U, - C/,+1 - (C/,_i - U,) = + 2U, - C/,_i, 

lo^^L- 2Uo, 

and 

j 

1=1 

= L-2ij + l)Uj+2jUj+i ^^°M^nj = 1). 
RecalUng M = —L/2, we obtain 

= (j = 1,2,..., 5-1). 

And, by definition of (^j f ([77|) ). we derive 

xp _ 6i _ L/q - 2^1 + j72 _ 2^1 2/1^0 



By poo]) , it follows that 
The fact that ([55]) ^s = iV = g + l completes the proof under the condition 

(IMl). 

For general cases, let us define g{i) (i = 1, 2, . . . , s) by 

= • • • = ^^(i) > '?e(i)+i = ■ • • = ^eia) > '^e(2)+i = • • • > <^e{s-i)-i = ■ ■■ = <ie{s)^ 
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and q{0) := 0. We obtain 
and 

Nj = lo + 2^ m {Li - ipQ+i) ) 
1=1 

s 

= + E (^>(fc)-^'(fc-l))(L^-L,(,+l)) 

fc=i e(fc-i)+i<i<e(fe) 

= L- 2{e{j) + !)[/,(,) + 2^.(J)c/,(,)+l. 

We can complete the proof in a similar manner to that of the previous case. I 

Acknowledgment. The authors are very grateful to Professor Ralph Willox for 
helpful comments on this paper. 

A The proofs of the remained lemmas 
A.l Proof of lemma [3.121 

To prove lemma [3.121 we investigate the 10-elimination (Section 12. 3p in detail. 
Let us consider a state consisting of N blocks of consecutive I's and O's, which 
are arranged alternatingly. We denote the length of the fc-th block of consecutive 
I's by Qk and fc-th consecutive O's by Wk- In order to show how these blocks 
are reconstructed by 10-eliminations, it is convenient to draw a graph which 
consists of nodes and links. For example, let us consider a state where = 3, 

2, Q3 = 6, and = 12. (See figure [H) The 

110000000 

T 

11100000000 

T 

110111000000000 

T 

1110011110000000000 

T 

11110001111100000000000 

\ 10-elimination 
1 1 1 1 10001001 1 1 1 1 1000000000000 

Figure 9: A state of pBBS with iV = 3 and the associated graph. The 
nodes at the bottom the graph are associated with (Qi, M^i, (32, M^2, Q3, W3) = 
(5,3,1,2,6,12). 

blocks of I's are represented by white nodes, and those of O's by black nodes. 
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A number is associated with each node. The numbers in the bottom of the 
graph are equal to Qi,Wi, . . . ,Qn, Wn respectively. Going up we arrange by 
one step, these numbers decrese by 1. The sign means zero, and a blocks 
of consecutive I's or O's disappears at the point where * appears. When one 
block disappears, the two blocks adjacent to a join together. Figure [TUl show 
examples of the graphs associated with a typical state. 

1111000000 

T 

1111100000000 

T 

1111110010001000000 

I 10-elimination 
11111 1 10001 100001 10000000 

11100 

T 

111101000 

T 

1111100110000 

T 

11111100011100000 

Figure 10: Examples of the graph associated with two typical states. In the 
first case, two blocks of O's disappear simultaneously. In the second case, two 
adjacent blocks disappear simultaneously, which would require writing the * 
twice. We shall write it only once. 

Let us define several terms relating to this associated graph. 
Definition A.l A tree is a connected component in the associated graph. 
Note that any tree has exactly one *. 

Remark A.l Only two types of tree can exist. One is a tree consisting of 
white nodes, and the other is a tree of black nodes. We denote the 'white tree ' 
as a'w-tree', and 'black tree' as a 'b-tree'. 

Definition A. 2 Let P be a node, and t be a tree in the associated graph. 

The height of P, denoted by Ilt(P), is the number of links in the path from P 

to the bottom of the graph. 

And the height oft, denoted by }lt(t), is the height of * contained within t. 

Let us denote by the graph associated with the state x G Z^. We 
introduce a semiordering on the set of trees by 

t,t' e <f>a;, t)~ t' ^ t straddles t' . 

Now we define two important sets. 
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Definition A. 3 Let t E be a tree, then we define a set of trees as 
Und(t) {s e ■ tree] s ^ t}, 

and a set of integers as 



is an associated number with 
a node at the foot of tree t. 



with X = {A^(i)}ffi, A21-1 = Qi, and A21 = Wi. 




h Wi Q2 W2 Q3W3 



t2 



Figure 11: An example clarifying Def lA.2l 

Example In figure llli ti and ^4 are white trees, and t2 and t^ are black trees. 
The trees have the relation ^4 ^ ^3 ^ ti, Und(ti) = {^3, ti}, Ft(ii) = {Qi, Q3}, 
etc. The height of ^3 and ^4 are given as Ht(t3) — 4, and Ht(t4) — 1. Note that 

t is a w-tree <^ Ft(t) C {Q^jili, (101) 

and 

t is a b-tree ^ Ft{t) C {W,}fLi- (102) 

Lemma A.l Let t & be a tree. The number of links in t is equal to 

A<,(i)GFt(t) 

Proof. This is a natural consequence of the definition of the associated graph. ■ 

Remember that a maximal element 1^ in a semiorderd set X is the element 
which satisfies 

r e A, ^ ^ e = e'- 



Remark A. 2 Let t Cz ^x be a w-tree, and s Cz ^x be a b-tree. Then any maximal 
element o/Und(t) \ {t} is a b-tree, and any maximal element o/Und(s) \ {t} is 
a w-tree. 
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We call the node which is connected with more than 2 links the branch point. 
We define the multiplicity of the branch point P as 



mp := {the number of links connected to _P} — 2. 

Note that 

{the number of maximal element in Und(i)\{i}} = mp. 

Pibranch pt. int 

Example In figure [TUl the multiplicity mp is equal to 2. The maximal elements 
of Und(i) \ {t} are t2 and tg. Note that Ht(P) = Ht(t2) = Ht(t3) = 2. 
Let us define 

Ft(Und(t)) := IJ Ft(s), 

sGUnd(t) 

where t £ \s a. tree. 

Lemma A. 2 Let t £ he a w-tree, and s £ <^x he a h-tree. Then, 



and 



t'eUnd(t) 

E Ht(s')^ 

s'GUnd(s) 



E 

Q,eFt(Und(t))n{Qj} 

E 

WieFt(\}nd{s))n{Wj} 



(103) 



(104) 



Proof. Let h := ^Undlt). We prove the lemma by induction of h. 

When ft, = 1, a w-tree t is a straight segment of length Qj, where {Qj} = 
Und(t). Then, Ht(t) ~ Qj. Hence, recalling (|10ip . we can conclude that (|103p 
is true immediately. Similarly, (I104p is true for a b-tree s. 

Suppose that (|103p and (I104p hold for h — 1,2, ■ ■ ■ ,p. Let i be a w-tree with 
T^Und(i) = p+1. Denote the branch points which belong to i by Pi, P2, • ■ ■ , -Pr- 



m{t) 
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Figure 12: An example of tree t. (7^ :— Ht(P,)). 
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From lemma lA. 11 we obtain 

{the number of links consisted in t} — Qi 

QieFt(t) 

r 

= Ht(0+^Ht(P,). (105) 

(See figure IT^) Since we can assume that the maximal elements si, S2, ■ ■ ■ , Sr G 
Und(t)\{t} satisfy Ht(sj) = Ht(Pj) (j = 1,2, . . . ,r), pOS)) becomes 

r 

0, = Ht(t)+^Ht(s,). (106) 
QieFt(t) j=i 

On the other hand, we easily find that 

r ^ I \ 

^Ht(.,)= mt')-Y.\ E mt")\. im 

J = l t'GUnd(t)\{t} J = l \t"eUnd(sj)\{sj} / 

Let tjk (fc = 1, 2, . . . , Ij) be maximal elements of Und(sj)\{sj}. By the induction 
hypothesis and (|103p we obtain 

E( E Ht(t")] = EE E Ht(t") (108) 

= EE E (109) 

j=l fc=l Q,eFt(Und(tj;.))n{Q„} 

E (110) 

Q.e(Ft(Und(t)\{t})) 

Substituting (fT07|) and pIO)) to (fT06l) . we find that p03)) holds for /i = p + 1. 
Equation (|104p can be proved in a similar manner. ■ 

Let 

Hk := {the height of the fc-th smallest tree in ^x}- 

Note that the proof of lemma 13.121 is completed by proving the following two 
formulae pTTjl and pT^ . 



mm ■ 



1 + ■ • • + fffe 



Y \=Hi + --- + Hk. (Ill) 

the number of boxes 

from under to the fc-th step in the Young diagram 

(112) 
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Proof of Let {Aa-(i)}* be an element of B{N — k, N) which satisfies 



mm 



{A<,(i)}ee(Ar-fc,Ar) J {A„(i)}*ei3(A'-fc,JV) 



Without loss of generality, we can assume Qa G {^^^(i)}* for some a. Let i £ 
be a tree with G Ft(i). It follows that 

G Ft(Und(t)) ^ G (113) 

("." Otherwise, there exists Q^jj , . . . , Q^, with Qp, G Ft(Und(i)) and Qp, ^ 
{Acr{i)}*. Let be the nearest branch point to Qp^- Without loss of gen- 
erality, we can assume Ht(Pi) = min (Ht(Pj)). Let us define the subgraph 

the path Pi ^ P consists of links 
extending below Pi 



<E>; := <! P : node 

and 



$** := {the node and links straddled by $*. }. 
Note that $** \ {g^J C {A„(^)}* ■ The inequality 

leads to a contradiction with the definition of {^^.(i)}*. ) 
The relation l|113p claims 

{A^i.)}* = |lJ(Ft(Und(t,)) n {Q„})| n |lJ(Ft(UndGs™)) n {Wn}) 
for some w-trees tj and b-trees Sm- From lemma [^31 we obtain 

{A„^i-^YeB{N-k,N) 

for some r. It completes the proof of (|llip . ■ 

Proof of (|112p Note that the number of boxes below the fc-th step in the Young 
diagram is equal to the number of bridges connected to the blocks which dis- 
appear by 10-eliminations, when the number of solitons becomes N — k. (See 
figure [5l) By definition of the associated graph ^x, it is clear that the equation 
(fn^ holds. ■ 

From proposition 13. 9[ (jllip and (|112p . Lemma [3. 121 is proved. 
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A. 2 Proof of proposition 13.71 and 13.101 

In this subsection, we give the proof of proposition 13.71 and 13.101 Recall : 

yri+l = (A - a„)j/„ - bnVn-l 

Let Xj ~ \ ~ aj and Yj ~ —bj. For convenience, we denote Xj by (j), and 
Yj by (j — in this section. For example, {a){b — l,b) means Xa ■ Y\j. 
Let 



and the initial condition 



{ii,i2, . . . ,is, ji,i2, ■ ■ -Jt} 



2s + t = k + l- j, 

{ii - . . . - i,«s,ii, ■ • ■ ,it} 
= {j,j + i,---,jfe} 



We define 

((j, j + 1, . . . , /c)) («1 - 1' ^l)(*2 - 1, *2) • ■ • («s - 1, ls){jl){j2) ■ ■ ■ (jt). 

{ii,i2,...,ia Ji,... Jt}GOjfc 

For example, ((1,2,3)) = (1)(2)(3) + (1, 2)(3) + (l)(2, 3) = XiX2X3+r2X3 + 
X1Y3. We also defined ((2, 3, . . . , k)) by the similar rule. (For example, ((2, 3)) = 
(2)(3) + (2,3) = X2X3 + r3.) 

Lemma A.3 If2<n<N + l, then a;„ = ((1,2, ...,n~ 1)). And if 3 < n < 
N +1, then y„ = Yi x ((2,3, . . . ,n - 1)). 

Proof. We prove the lemma by induction. Whenn = 2, a;2 = (A— 01)2:1— 61X0 = 
A — ai = Xi ~ ((1)). Since Xn+i = XnXn + YnXn-i, from the induction 
hypothesis, x^+i = (n) • ((1, 2, . . . , n - 1)) + (n - 1, n) • ((1, 2, . . . , 71 - 2)) = 
((1, 2, . . . , n)). Hence Xn — ((1, 2, . . . , n)) holds for any n. We can prove the 
assertion for y„ in a similar fashion. ■ 

If ft > iV, n should be considered as an element of Z/iVZ because of the 
periodic boundary condition a^+j = aj, b^-^-j = bj. In this sense, we may 
write Yi = (A^, 1). Propositions 13.71 and 13.1171 are proved immediately from 
lemma 1X73] In fact, the relation yN+i = Yi x ((2,3, . . . ,N)) is equivalent to 
proposition 13.101 In order to prove proposition 13.71 note that A(A) = xn+i + 
VN = ((1,2, ...,Ar)) + (Ar, 1)((2,3, ...,7V-1)). Let us decompose ^(7V) C 2^/^^ 
into two disjoint sets 

A{N) ^Mu A' 

where consists of the subsets (c 7j/N7L) that include (A^, 1), and A" := 
A — A' . Proposition 13 . 71 is obtained by the fact that 

^ Y,, . . . Y,,X,, . . . X,„_,, = (iV, 1)((2, 3, . . . , TV - 1)) 
(ii-i ji,--- jfc-i 
and 

^ y,,...y,,x,,...x,„_,, = ((i,2,...,iv)). 

(ji-i ji>--- jfc-i 
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A. 3 Calculation of Ej 

Recall that aj (j = 0, 1, . . . , g) are the roots of A(A) — 0, and {j — 1,2, ... , g) 
are the roots of j/Ar+i(A) = 0. Hence, 

A, A(a^,) _ m,2,...,N)) + {N,l){{2,Z,...,N~l))]^=^^ 

where iV = g + 1 . 
Since 

((2,3,...,iV))U=^, =0 

and 

((1, 2,...,N)) = (1)((2, 3, . . . , iV)) + (1, 2)((3, A,..., N)), 

we find 



YT, , [(1, 2)((3, 4, . . . , TV)) + {N, 1)((2, 3, . . . , iV - 1))]a. 

Note that (j) = A - + Fj) and {j - = -/j-Vj. 

We show some examples here. 
Example 
WW~2{g^l), 

loVo + hVi u loVo + hVi 



-(^1 - ai) = 



Ml - "0 Ml 

^> - min {Qj] + + - (-Pi) 



If iV = 3(.g = 2), 

n (M. ~ g.) = ' + ' ^^^^ ' + "^^^^ . (114) 

To calculate the ultradiscrete limit of the right hand side of (|114p . we have to 
compare the magnitude of each term. 

For example, when j = 1, (P2 - A) < min {Q\, W^), and Ql + + P2 - A < 
Q\ + + min {P2 - Pi, Qg, M^2°}, then 

Si =Q!] + W^o + P2-Pi. 
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